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This article extends previous studies on the effect of tensor terms in the Skyrme energy density 
functional by breaking of time-reversal invariance. We have systematically probed the impact of 
tensor terms on properties of superdeformed rotational bands calculated within the cranked Hartree- 
Fock-Bogoliubov approach for different parameterizations covering a wide range of values for the 
isoscalar and isovector tensor coupling constants. We analyze in detail the contribution of the tensor 
terms to the energies and dynamical moments of inertia and study their impact on quasi-particle 
spectra. Special attention is devoted to the time-odd tensor terms, the effect of variations of their 
coupling constants and finite-size instabilities. 
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I. INTRODUCTION 

Recent years have seen a renewed interest in the role 
of the effective nucleon-nucleon tensor force for nuclear 
structure, sparked by the finding that it provides one 
of the possible sources for the evolution of nuclear shell 
structure with neutron and proton numbers. Indeed, the 
contribution of tensor interactions to single-particle ener- 
gies depends on the filling of shells. It (nearly) vanishes 
in spin-saturated nuclei, whereas it might be significant 
when only one out of two spin-orbit partner levels is filled 
for one or even both nucleon species [T]. 

Up to now, none of the standard parameterizations 
of any mean-field approach considered an explicit ten- 
sor part, cf. Ref. |2] for a historical overview. The first 
studies of the effective tensor interaction within self- 
consistent mean-field approaches concentrated on single- 
particle spectra in chains of semi-magic spherical nuclei 
covering all successful models, i.e. the non-relativistic 
Gogny force [3l S] and Skyrme interactions [21 UHl] , as 
well as relativistic mean- field approaches [51 [S] ■ More re- 
cently, the impact of the tensor terms on more complex 
structure properties has been studied as well, such as 
the topography of deformation energy surfaces [10] and 
various spin- and spin-isospin excitation modes in Quasi- 
particle Random Phase Approximation (QRPA) using 
Skyrme functionals [111117] and Gogny interactions [4 . 

These QRPA calculations deal with a very different 
aspect of an effective tensor interaction than the analy- 
sis of single-particle energies. This becomes most obvious 
when using the Skyrme energy density functional (EDF). 
The Skyrme EDF can be separated into two parts: the 
first one composed of densities and currents that are even 
under time reversal such as the normal and kinetic den- 
sities, and a second one grouping combinations of time- 
odd densities such as spin density or current. The latter 



part of the EDF is usually called the "time-odd" part, 
although strictly speaking the EDF itself is time-even by 
construction. While these time-odd densities are zero for 
the HFB ground states of even-even nuclei, they become 
non-zero for 

(i) blocked quasiparticle states, i.e. the self-consistent 
calculation of non-collective low-lying states in odd- 
A and odd-odd nuclei, or K isomers in even-even 
ones. The time-odd terms contribute to the total 
energy [7] and their presence can strongly 
modify the expectation values of time-odd observ- 
ables such as magnetic moments [Ml ES] , 

(ii) rotational states calculated by the cranked HFB 
method [25115^ , where they affect the alignment of 
single-particle levels with the rotational axis and 
thereby the moments of inertia, 

(iii) time-dependent Hartree-Fock (-Bogoliubov) 
(TDHF(B)) [33H35] and its linear response hmit, 
the (Quasiparticle) Random Phase Approximation 
((Q)RPA), 

(iv) configuration mixing such as symmetry restora- 
tion or Generator Coordinate Method calcula- 
tions, or adiabatic Time-dependent Hartree-Fock- 
Bogofiubov calculations (ATDHFB) [551|5S1I57]. 

The study of excitation modes of unnatural parity, such 
as for example Ml, spin-dipole or Gamov- Teller exci- 
tations in (Q)RPA provides a sensitive benchmark for 
the time-odd terms in the EDF. Indeed, for those the 
residual interaction is entirely determined by the time- 
odd terms, cf., for example, Refs. [551H5] and references 
therein. For the other phenomena listed above, the time- 
odd terms provide a correction to the dominant time-even 
terms that might become substantial in some cases. One 
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such observable are the moments of inertia at high spin 
in superdeformed (SD) rotational bands of heavy nuclei 
[33H1SI- These will be the object of our study. 

Skyrme's two-body tensor force contributes to the 
time-even and time-odd parts of the EDF. Studies of the 
eigenvalue spectrum of the single-particle Hamiltonian of 
even-even nuclei only probe the contribution to the time- 
even part. The corresponding time-odd terms affect how 
the nucleus responds to its collective rotation, thereby 
modifying its moment of inertia and how it evolves with 
spin. Of course, the time-even tensor terms also influence 
the moments of inertia. 

The aim of the present study is to investigate the 
generic influence of tensor terms on high-spin properties. 
The following questions will be addressed: 

• How do the time-odd tensor terms behave when 
increasing the total spin of the nucleus? 

• How does the presence of time-even and time-odd 
tensor terms influence the dynamical moments of 
inertia in superdeformed rotational bands at high 
spin? 

• How much of these changes is caused by the time- 
even part of the EDF, i.e. the modiflcation of 
the single-particle spectrum at spin zero, and how 
much by the time-odd part of the EDF? 

• How much of these changes is caused by the tensor 
terms themselves, and how much is caused by the 
rearrangement of all other terms during the fit of 
the parameterizations? 

Studies of the impact of time-odd terms on the moments 
of inertia in superdeformed bands have been performed 
before in the context of Skyrme interactions [2Sj |40j and 
relativistic mean-field Lagrangians [301 133 , but none of 
these studies considered time-odd terms associated with 
genuine tensor interactions. 

The present article complements the studies of spheri- 
cal single-particle energies of Ref. [5] and of the deforma- 
tion energy curves of Ref. (TUI. We will refer to these 
references as Articles I and H. The present Article is 
structured as follows: In Section [H] we briefiy review 
the properties relevant for our discussion of the Skyrme 
EDF including tensor terms. In Sect. |III[ we analyze 
in detail how sensitive is the superdeformed rotational 
band in ^^''Hg when the coupling constants of the Skyrme 



EDF are varied. In Sect. |IV[ we check the generality 
of our conclusions for ^^"^Hg by similar calculations for 
the SD ground-state band of ^'"^^Dy, and Sect. |v] sum- 
marizes our results. Appendices provide further techni- 
cal information about the interrelations between the cou- 
pling constants of the Skyrme EDF (Appendix |A]), the 
representation of local densities and currents in our code 
(Appendix [b]) and about the Landau-Migdal interaction 
corresponding to the standard Skyrme EDF with tensor 
terms (Appendix ICl) . 



II. THE SELF-CONSISTENT MEAN-FIELD 
METHOD 

The energy of the atomic nucleus can be expressed by 
means of an energy density functional |47H50) , containing 
five parts: the kinetic energy, a Skyrme potential energy 
functional modeling the strong force in the particle-hole 
channel, a pairing energy functional, a Coulomb energy 
functional, and terms to approximately correct for the 
spurious-motion caused by broken symmetries 



S — £ kin + ^Sk + ^1 



-'pairing 



Coulomb 



(1) 



For the kinetic energy and the Coulomb energy func- 
tional comprising a direct term and the exchange term 
in Slater approximation we use the same expressions 
as presented in Ref. |35]. For all parameterizations 
used throughout this Article, the center-of-mass recoil 
effect is approximately taken into account by subtract- 
ing £corr — {TlikPk) /'^''^^ from the total energy, which 
amounts to an A-dependent renormalization of the nu- 
cleon mass. In the following Sections we will introduce 
the ingredients of the Skyrme EDF, its explicit form and 
the equations of motion. 



A. Densities and currents 

Under the assumption that the single-particle states 
are either neutron or proton states, the Skyrme part of 
the energy density functional £sk depends on the follow- 
ing local densities and currents 
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which are the density Pq{r), the kinetic density Tq{r), the 
spin-current (pseudo-tensor) density Jq,^,y(r), the cur- 
rent (vector) density jg(r), the spin (pseudo- vector) den- 
sity s^(r), the spin- kinetic (pseudo- vector) density Tq(r), 
and the tensor-kinetic (pseudo- vector) density Fg(r) for 
protons and neutrons {q — n,p). All densities and cur- 
rents can be recoupled to isoscalar t — and isovector 
(t = 1, is = 0) densities (e.g. po{r) = pp{r) + p„(r) and 
Pi{r) — Pn{r) — Pp{r)). They are constructed from the 
density matrix in coordinate space and its derivatives up 
to second order [501 EH 

Pq{r,a,r'a') = (a^(r', cr', g)a(r, ct, g)) , 
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ipg(r, r')5aa' + ^{a\6-\a') ■ Sg(r, r') 



with 



(3) 



(4) 



here to aspects relevant for this work. As outUned in the 
Introduction, the Skyrme EDF can be separated into two 
parts 



a=±l 



" + n 1 



(12) 



t=o,i 



89(1", I"') = X! Pgi^'^'^''^') i'^'l^W) ' (5) a part that contains only time-even densities 



(Ta-'=±1 



and (T are the Pauli spin matrices. From the properties 
of the density and the spin density matrices under time- 
reversal EH 



s^r,!-') = -s,(r',r) 



(6) 
(7) 



it follows that 



s[(r) - -s,(r) , j^(r) = -j,(r) , T^(r) = -T,(r) , 



F^(r) = -F,(r) 



(8) 



One notes that Pq{r), Tg(r), and Jq,pi/(r) are time-even, 
whereas Sq(r), jg(r), Tq(r), and Fq{r) are time-odd. 
When time reversal is a self-consistent symmetry, the 
time-odd densities vanish. When intrinsic time-reversal 
invariance is broken, both time-even and time-odd den- 
sities contribute to the energy and the single-particle 
Hamiltonian. 

For the part of the EDF that describes the pairing 
correlations, ^pairing, we also need to introduce the skew- 
symmetric pair tensor [52 

Kq{r, cr, r', a') = (a(r', a', q)a{r, a, q)) , (9) 

or, equivalently, the pair density matrix |53j 

Pg{r, cr, r', cr') = -cr'(a(r', -cr', q)a{r, cr, q)) 

= -a'K,(r,a,r',-(T'). (10) 

The pair density matrix presents the interest that it en- 
ables to construct a local pair density 



Pg(r) 



cr = ±l 



Pq(r,cr,r,cr) , 



(11) 



that facilitates the construction of a local pairing EDF 
as used here. The pair density pq{r) is neither time-even 
nor time-odd and becomes complex when intrinsic time- 
reversal symmetry is broken. 



B. The energy density functional 

1. The Skyrme energy functional 

Skyrme functional have been discussed extensively in 
the literature Pl [?n 155 1 liMCT l [5i] . We restrict ourselves 



CMpi + Ct'pt^pt + Clptn + G^ -'ptV ■ J 



(13a) 



of time-odd densities 

nl °- = QiPoH + c7^(v • + c^^s, ■ As, ~ cu't 

+ C^st ■ Tt + Cfst ■ Ft + CY 'st • V X ji . 

(13b) 

The common practice to call these energy densities 
"time-even" and "time-odd" refers to the properties un- 
der time reversal of the densities and currents they are 
built from and not to the properties of the energy density 
itself, as the EDF is time-reversal invariant by construc- 
tion l55l. 



All coupling constants of ( 13a I and ( 13b ) could be cho- 
sen to be density dependent. ^ In practice, however, the 
density-dependence is usually restricted to the Cf and 
coupling constants, and chosen to be a non-integer 
power of the isoscalar density [SH [S5] 



cM^cm + {c?[pnn.]-cm 
cnpo] = cm + (cnpnn.]-cm 



Po 
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Po 
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(14) 
(15) 



where pnm is the value of the isoscalar density po in sat- 
urated infinite nuclear matter. 



2. Choice of independent coupling constants m the EDF 

The Skyrme energy density functional (EDF) can 
be introduced in two non-equivalent ways (49) . One 
can start from a density-dependent zero-range two-body 



^ If all coupling constants are density-dependent, there appear 
additional terms containing the nabla or Laplacian operator act- 
ing on a local density, as some bilinear forms are not equivalent 
under partial integration anymore |55| . 
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interaction as proposed by Skyrme [571 ISH] and de- 
fine the Skyrme EDF as its Hartree-Fock expectation 
value (231 [Sni IMl- We refer to this approach as force- 
generated EDF. It is customary that the pairing func- 
tional is calculated from a different pairing interaction, 
such that only the direct and exchange terms in the 
particle-hole channel refer to the same interaction, but 
not those in the pairing channel. The same expression for 
the energy densities T-Ll '^' and Hl °' is obtained when writ- 
ing down all possible bilinear combinations of the local 
densities (2a)-(2g| up to second order in the derivatives 



that are invariant under parity, time reversal, rotations 
and Galilean transformations [53 ES]. Then, the coupling 
constants of most time-odd terms can be chosen indepen- 
dently of those of the time-even terms. By contrast, in a 
force-generated EDF, the coupling constants of time-odd 
terms depend on those of the time-even ones; the 9 free 
parameters of the most general density-independent two- 
body Skyrme force determines all 18 coupling constants 
C of the corresponding EDF ([l3a][l3b|. The various pos- 
sible choices for the coupling constants of time-odd terms 
that correspond to the same time-even part of the EDF 
have been reviewed in Article II. 



Galilean invariance is a necessary constraint to obtain 
a functional that depends only on the relative momenta 
of the nucleons, but not on the total momentum. This 
invariance is particularly important for dynamical calcu- 
lations such as TDHF and TDHFB [31] or cranked HFB 
such as performed here [5M5T] since it ensures that the 
results of the calculation will not depend on the frame 
of reference. In relativistic approaches, the EDF has to 
be constructed to be Lorentz invariant instead [30, 32]. 
Galilean invariance is automatically fulfilled for an EDF 
generated from the Skyrme force. In a functional ap- 
proach, it is automatically fulfilled for some terms and 
has to be imposed by taking specific combinations of 
other terms [2HI 1331 • This results in the presence of 



3. Skyrme's tensor interaction 

A zero-range tensor force with two terms was originally 
proposed by Skyrme [57l [58] ^ 



v\r) = |{ [3(^1 • k')(<T2 • k') - [a, ■ a,)k'']5{r) 

+ S{r) [3{CTI ■ k){cT2 ■ k) - (^1 • <T2)k2] } 

+ |{3(<Ti • k')<5(r)(<T2 • k) - (^1 • ^2)k' • ,5(r)k 

+ 3(*T2 • k')<5(r)(cri • k) - (o-i • ^2)k • J(r)k'} . 

(16) 

The inclusion of this tensor force gives rise to new 
terms in the force-generated energy density [2] 

z 



C 



St 



Jt., 



fl—X 



2 5^ Jt,^yJt,u^ 

(17) 



fl^V — X 



The first two terms also appear in the EDF constructed 
from a central Skyrme force, whereas the latter two occur 
for genuine two-body tensor forces only. They differ in 
the way derivatives and Pauli matrices are coupled. In 
the first two terms, the scalar products are between one 
derivative and the other, and between one Pauli matrix 
and the other. In the last two terms, the scalar products 
are between derivatives and Pauli matrices. 

In a force-generated EDF, each coupling constant 
and C^" results from two contributions. Following Ar- 
ticle I, we label the contributions coming from the cen- 
tral Skyrme interaction by the letter A and those gen- 
erated by the tensor part of the interaction by B. The 
inclusion of a tensor force thus increases the flexibility 



CI, CY'-', , and Cf in both the time-even and the for the choice of the coupling constants Cj' and Cf" in 



time-odd energy density ( 13a|13b ). Since several time- 
odd terms of an EDF are not constrained by symmetry 
requirements, the number of independent coupling con- 
stants is always smaller for a force-generated EDF. This 
renders the mapping of an arbitrary local EDF onto a 
density-dependent Skyrme force a priori impossible when 
time-reversal invariance is broken. 



a force-based approach. At the same time, the tensor 
force introduces additional terms in the EDF that couple 
derivatives and Pauli-spin matrices in a unique manner. 

The most appropriate way would be to label tensor 
terms as those generated by the tensor force, Eq. (16 1. 



However, the above discussion shows that these terms 
cannot be easily singled out with respect to similar terms 
generated by the central part of the EDF. Therefore, 
throughout this article, we will call "tensor terms" those 
terms in the EDF that couple two Pauli matrices and 



In the literature, there also exist hybrid approaches 
which adopt the larger freedom of functionals for some 
terms in the functional only and follow the "force- 
generated" philosophy for all others. We will come back 
to this in Sect. IIIFI 



^ Skyrme's tensor force is sometimes given in a different, slightly 
simpler form [2l 1101 162H64I . Both forms give rise to the same 
EDF, but not to the same residual interaction, for example, in 
QRPA. 



5 



two derivatives. Although not all of them are related to 
a two-body tensor force, they are of the same order as 
the terms bilinear in the spin-current tensor density that 
have been called tensor terms in Articles I and II. 



4. The pairing energy 

As in our previous studies, we have chosen a density- 
dependent zero-range interaction to describe the pairing 
correlations |49l ESI [65] , which leads to a functional of 
the form 



pairing 



4 

q—p,n 



Po(r) 



p,(r)p;(r) , (18) 



where the switching density pc determines whether the 
pairing is more active in the volume of the nucleus or 
on its surface. The functional depends on the local pair 
density pq and the local density pq. 



C. The cranked HFB method 

The superdeformed rotational bands will be calcu- 
lated by the self-consistent cranked HFB approach. This 
method can be seen as a semi-classical description of the 
collective rotation of a finite system with a constant an- 
gular velocity oj. In particular, it takes into account the 
distortion of the nucleus' intrinsic state due to the cen- 
trifugal and Coriolis forces that are induced by the col- 
lective rotation [5^ [551 IMl . 

The variation of the EDF including constraints on par- 
ticle number, orthonormality of the quasiparticle states 
and rotational frequency leads to the cranked HFB equa- 
tion 



h — X — io^J^ 
-A* 



A 

-I- A + uJzJt 



E,. 



(19) 

where C/^ and are the two components of the quasi- 
particle wave functions and the quasiparticle energies, 
often called Routhian in this context. The effective in- 
teraction enters the HFB Hamiltonian through the mean- 
field Hamiltonian h and the pairing field A 



~ Sp- 



(20) 



The Fermi energies A for protons and neutrons and the 
rotational frequency in Eq. ( 19 ) are the Lagrange mul- 



tipliers of the constraints, which are self-consistently ad- 
justed to fulfill auxiliary conditions for the mean values 
of the particle number and of the projection along z of 
the angular momentum. The component Jz of the angu- 
lar momentum J is chosen along the axis perpendicular 
to the axes of longest elongation. At high spins and large 
deformation, the solution of Eq. ( 19 ) can be shown to be 



an approximation of a variation after projection on an- 
gular momentum |68) . As such, the model is particularly 
well adapted for the description of superdeformed bands. 

All calculations have been carried out using the tri- 
axial self-consistent cranking code CR8 documented in 
Refs. [271 [Ml [69]. The HFB equation is complemented 
by the Lipkin-Nogami (LN) prescription to avoid a sud- 
den breakdown of pairing correlations as a function of 
rotational frequency. 

We recall that for constrained calculations, as dis- 
cussed below, the constraints do not contribute to the 
observable total energy, whereas the eigenvalues E^^ of 
the HFB Hamiltonian used to construct the quasiparti- 
cle Routhians contain a contribution from the constraint. 



D. The single- particle Hamiltonian 

The isospin representation of the EDF is the most ap- 
propriate one to discuss its physics content. However, 
a representation where proton and neutron densities are 
explicitly used is more convenient for numerical imple- 
mentations. In this case, one possibility to write the 
Skyrme EDF is [27 



H{v) -t- Uqiv) 



(21) 



q—p,n 
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n{r) = 61 p2 + 63 (pr - f) + fes pAp + 67 
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s T 



z 2 z 
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(22) 



containing the total density p = Pp + Pn and similar for 
all the other local densities and currents, ^ and 

Hqir) ^b2pl+ bi {pqTq - j^) + be PqApq + feg pQ P^ 

+ bgq {pqV ' Jg + j,, ' V X S,) + 611 + &13 p"s^ 



z 



^ ^ Jq^fiuJq^u^ 2 Sq • Fg 



^ Even though the "total" local densities and currents are identical 
to the "isoscalar" local densities and currents, we use a different 
notation to clearly distinguish between the isospin representation 
and the proton-neutron representation used in our codes. 
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+ big Sg ■ ASq + 621 (V • S,)' 



(23) 



containing the proton and neutron local densities and 
currents. The relation between the coupling constants 



in the isospin representation, Eq. (13a I and (13b), and 



the parameters in the proton-neutron representation, 
Eq. (21), is given in Appendix [Aj The single-particle 
Haniiltonian for protons and neutrons is given by 

hg^-V- Bg{r)V + Ug{r) + Sg(r) ■ & 

z 

-^[A,(r).V + V.A,(r)] 

- V- [(T-C,(r)]V- V-D,(r)o-- V. (24) 

The expressions obtained in [57j for the inverse effective 
mass Bg{r), the single-particle potential Uq{r), and the 
time-odd field Ag(r) are not affected by the introduc- 
tion of tensor terms. The local potentials that contain 
contributions from the tensor terms are given by 



y , EkAII/ (69 + bgg V^Pg) 

K — X 

2 614 J^j^ H- 2 615 Jq^^u 
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2 61 



K — X 

z 

K — X 



(25a) 



(69 V X j + 695 V X j 



+ 2610 +2611 Sg^fj, 

+ 2 6i2/5".s^ + 2 6i3 Sg^^ 

— bi4 — 615 Tg^fj^ 

— 2 616 — 2 617 Fq^^ 

+ 2bisASf, + 2biQAsg^f, 

— 2 620 V^(V • s) - 2 621 V^(V • s,) , (25b) 
Cq,p(r) = - 614 - 615 s,,^ , (25c) 
Dq,f,{r) =-2bieSf,-2bnSg^i,, (25d) 

The scalar central potential Ug{r) = S£/Spg{r), 
the position-dependent inverse effective mass Bg{r) — 
h'^/2m*{r) = S£/6Tq{r), and the spin-current tensor po- 
tential Wg^f^^{r) = 5S / SJq fj^i/ir) are all time-even fields, 
whereas Ag(r) = ^f/Jjg^r) and Sq(r) = 5£/5sg{Y), 
Cg{Y) = 5£/6Tg{Y), and D<,(r) = 5£/6¥g{T) are time- 
odd fields. The vector potentials A^(r) and Sg(r) are 
nuclear counterparts of electromagnetic potentials that 
couple orbital movement and spin to magnetic fields. The 
field Cg(r) introduces a spin dependence of the position- 
dependent effective masses of protons and neutrons. Fi- 
nally, the field Dq(r) contributes to a non-diagonal tensor 
effective mass that is position- and spin-dependent. As 



long as time-reversal invariance is not broken, the time- 
odd fields remain zero. 

For density-independent , C]^ , and Cf , in a static 
calculation, and for our choice of symmetries (see Ap- 
pendix |b] on conserved symmetries) the single-particle 
Hamiltonian ( 24 ) can be reduced to 



hg=-V- Bg{v)V + Ug{r) + Sgiv) ■ & 

z 

-« W^,,p^(r)V^a^-iAg(r)- V 

- V- [6--Cg(r)]V- V-Dg(r)^- V. (26) 

As already mentioned in the Introduction, the con- 
tribution of the tensor interaction to the eigenvalues of 
the single-particle Hamiltonian depends on the filling of 
shells. This is due to the near cancelation of the contri- 
butions of two spherical spin-orbit partners to the spin- 
current tensor J^i, when both levels are filled. As a con- 
sequence, the time-even terms bilinear in J^^, (nearly) 
vanish in spin-saturated nuclei, whereas they might be 
quite large when only one level out of two spin-orbit part- 
ners is filled, cf. the discussion in Articles I and II and 
references given therein. It is noteworthy that none of 
the various time-odd terms associated with the tensor 
force has such property. Apart from the usual cancela- 
tion of Kramers-degenerate levels that are connected by 
time reversal, there is no additional dependence of these 
time-odd terms on shell structure as such. 

Further technical information about the detailed form 
and symmetries of the local densities and fields as imple- 
mented in the CR8 code is presented in Appendix [B| 



E. Landau parameters 

The gross properties of the spin-spin interaction in 
nuclei are often characterized by the so-called Landau- 
Migdal parameters HHMl gDl [701 [TT] • In Landau theory 
for normal Fermi liquids |72| , the Landau-Migdal param- 
eters represent the strength of the residual interaction 
between particles on the Fermi surface. Being a sim- 
ple number in each partial wave and spin-isospin channel 
of the central and tensor interaction, they cannot rep- 
resent all details of the effective interaction in a finite 
system. Still, they provide an often useful first indica- 
tion about its relative strength. Relevant for our present 
study are the Landau parameters in the spin- and spin- 
isospin channels of the central residual interaction and in 
the tensor channel 
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FIG. 1: Coupling constants Co and Cf , where Cf = — Ct^ + 
^Cf , for the parameterizations discussed in the article. 



(27) 



We use the convention of Refs. 5391 [40] where the normal- 
ization factor is defined as the average level density A'o = 
2kpmQ/h'^n^ at the Fermi momentum kp — (|7r^po)^^'^, 
in which TOq is the isoscalar effective mass associated with 
a given parameterization, but other choices are some- 
times found in the literature. All higher Landau param- 
eters are zero by construction for a Skyrme EDF that 
contains only terms up to second order in derivatives. 

In a force-based framework, the central and tensor 
parts remain separated in the residual Landau interac- 
tion, such that ho and /ig are entirely determined by 
and to of Eq. (161, which at the same time do not con- 
tribute to fi, gi or cf. the expressions given in 
Ref. 14 . In a functional-based framework as assumed 



in Eq. (27), however, this clear separation is lost, as can 



be seen from the appearance of in all six Landau 
parameters. The reason is that one has to combine con- 
tributions from the J^^J^^v, JtiuJ^fj., Jtiti, s • F, and s • T 
terms to recover the structure of the tensor operator that 
multiplies ho and /iq, cf. Appendix [C] for details of the 
derivation. 



F. Parameterizations 

In Article I, a set of 36 parameterizations for the 
Skyrme interaction including a zero-range tensor force 
has been determined using a fitting protocol almost iden- 
tical to the one used for the SLya; parameterizations 
[371 13H]- These parameterizations, labeled T/J, system- 
atically cover a wide range of the C/ = —C^ + \ cou- 
pling constants of the tensor terms in spherical symmetry, 
see Articles I and II for further details. In the present 
study, we restrict ourselves to a subset of four of these 
parameterizations, T22, T26, T44, and T62 (see Figs, [l] 



and [2]) . T22 has been constructed to give vanishing con- 
tributions of the tensor terms in spherical symmetry and 
time reversal invariance. It is aimed to have properties 
close to those of SLy4, which does not include a tensor 
interaction and for which the contributions of the cen- 
tral part of the interaction to tensor terms have been ne- 
glected. However, the tensor terms of T22 can be differ- 
ent from zero when spherical symmetry or time-reversal 
invariance are broken, cf. Article II for the breaking of 
spherical symmetry. In the same way, time-odd terms 
C]^St ■ Tj and C^St • do a priori not cancel each other. 

The T22 and T44 parameterizations have an isovector 
C( coupling constant equal to 0. The isoscalar coupling 
constant Cq has the same value 120 MeV fm^ for T26, 
T44, and T62, while the isovector G( coupling constants 
respectively take the values 120 MeV fm^, MeV fm^, 
and —120 MeV fm^. As a consequence, the tensor terms 
in spherical symmetry are purely between particles of 
same isospin for T26, purely proton-neutron for T62 and 
a mixture of both for T44. 



As will be discussed in Sect. IIIB the force-generated 
values of at least one of the Cf^md CY' coupling con- 
stants of all T/J parameterizations leads to an unphys- 
ical solution of the cranked HFB equation at high spin. 
Unless noted otherwise, we have set C^' and Cj" to zero 
for those parameterizations. 

As a reference without tensor terms, we have also per- 
formed calculations with the SLy4 parameterization [48] 
for which C]^ is put to zero during the fit. To fix the 
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time-odd terms, we have adopted a hybrid point of view 
that has already often been used in the hterature before 
[?7l [75H75] : a force-generated point of view is taken for 
the coupHng constant C*, whereas we use the functional 
point of view to set the coupling constant of the C^^'^SfAs 
to zero such that all three terms of second order in the 



derivatives and Pauli matrices vanish C^^ 



CJ = 0. 



This choice is not unique. Other groups use SLy4 in 
their cranked HFB or QRPA calculations by setting just 
= 0, keeping and Cf^^ at their Skyrme-force val- 
ues [niiHiiin]- 

For the the pairing EDF (18), we choose a surface- 
type interaction with pc — 0.16 fm'^ and a strength of 
Vq — —1250 MeV fm~'^ for both protons and neutrons, 
together with a 5 MeV cut-off above and below the Fermi 
level as explained in [75] . 



III. RESULTS FOR SUPERDEFORMED BANDS 
IN i"HG 

A. General comments 

In contrast to superdeformed bands in nuclei around 
A = 150, the high-spin properties of nuclei in the Hg 
region are sensitive to pairing correlations. In cranked 
Woods-Saxon and Nilsson-model calculations [78] . 
the gradual increase of the dynamical moments of inertia 
j/^^) as a function of the rotational frequency results from 
both the alignment of the intruder orbitals and from a 
gradual disappearance of pairing correlations. The prop- 
erties of the ground superdeformed (SD) bands have been 
studied extensively within self-consistent cranked HFB 
models using an effective EDF. In general, a very good 
agreement with experiment is obtained for the Hg region 
[551 IMl [75H5T] . In view of this success, we choose the 
ground SD band in ^^"'Hg as a laboratory for the study 
of tensor terms on high-spin properties. 

Our discussions are mainly based on the behavior of 
the dynamical moment of inertia JT^^^-* as a function of 
huj 



Ji2) 



doj 



(28) 



where (Jz) is the average value of the projection of the 
angular momentum on the rotation axis. The relevance of 
this quantity for the purpose of our study becomes clearer 
when realizing that is proportional to the derivative 
of the EDF with respect to rotational frequency [H7| [12 



J 



(2) 



CO duj 



(29) 



which allows to calcu late the contribution of each term 
in the EDF ( |l3a|l3b| to J^^^ separately. The moment 
of inertia can also be decomposed into the neutron and 
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FIG. 4: (color online) (a) Dependence of the Cq^^sq • Aso 
term of a variant of the T22 parameterization on the value of 
Cq^" for the Jz — 54h state in the ground-state superdeformed 
band of ^^"^Hg (see text), (b) Dependence of all other time- 
odd terms containing the spin density St relative to their value 



at C(\ 



in the same calculations. 



proton contributions to Jz using Eq. ( 28 1 . These con 



tributions, however, do not correspond to the decompo- 
sition of the EDF into neutron-neutron, proton-proton 
and proton-neutron terms. 

The numerical determination of JT^^^-* is far from be- 
ing trivial. To obtain a smooth dependence of J^^^ as a 
function of fku, it requires a very high degree of conver- 
gence of the calculations. The derivatives with respect 
to uj are determined by finite differences formulas. It 



should be noted that Eqns. (28) and (29) might lead to 
slightly different j/^^^ , especially when convergence to a 
very high degree is difficult to attain. In general, the 
( Jz) are converged to a higher degree than the energy £, 
hence Eq. (28) is expected to be more stable. 



The experimental value of hu is given by and 
the one of the dynamical moment of inertia by JT'-^-' = 
Afi^ /AEj, where AE^ is the difference between two suc- 
cessive 7-ray energies populating and depopulating a 
level. Note that both can be determined without an an- 
gular momentum assignment of the level, as is the case 
for the ground-state SD band in ^^"^Hg. 
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B. Finite-size instabilities 

In our calculations of the SD rotational bands with 
T/J parameterizations mentioned before, we systemati- 
cally encountered non-convergence of the code when us- 
ing force-based coupling constants for all time-odd terms. 
After careful analysis, this behavior turned out not to be 
a numerical problem, but a property of these parameter 
sets. Switching to a functional framework, systematic 
variation of coupling constants reveals that large posi- 
tive or negative values of either C/^" or C^^ lead to an 
unphysical finite-size instability of a given parameteriza- 
tion of the interaction. 

As a typical example. Fig. |4] presents the energies of 
time-odd terms containing the spin density for the SD 
Jz = 54h state in ^^^Hg as a function of the value of 
C^** for a variant of the T22 parameterization. Panel 
(a) displays the absolute energy of the C^* Sg • Asq term, 
whereas panel (b) presents the evolution of other time- 
odd terms that contain the spin density relative to their 
value at C^^ = 0. In this calculation, the coupling con- 
stants of all time-even and time-odd terms are set to their 
force-based values, except for C^^" = Cq'^ — CY" = 
which are set to zero, and C^^ that is systematically 
varied. For larger or smaller values of than those 
shown in Fig. [4| our calculations do not converge. When 
approaching C^* « 36 MeV fm'^, the energy of the 
C^^ So • Aso term displays a steep downwards slope. Si- 
multaneously, all other terms containing the spin density 
are strongly amplified, in particular the Cq Sg term, indi- 
cating a strong change in spin polarization. Still, in spite 
of their strong variation, the absolute contribution of all 
these terms to the total energy of ^^"^Hg remains less than 
0.3 % even at the threshold of the finite-size instability. 
Also, a strong dependence of the spin terms on the cou- 
pling constant is not a necessary condition for the onset 
of an instability of the C^^ Sq-Asq term. This can be seen 
when approaching C^^ sa — 24 MeV fm^, beyond which 
the interaction also becomes unstable. Similar results are 
obtained for the variation of the C^'*; an instability sets 
in at the same values as for C^^. In fact, the instabil- 
ity of the Cf^'^ St ■ Ast terms at large positive values of 
C^'^ ^ 36 MeV fm^ has already been pointed out earlier 

[211 IM]- 

In a similar manner, we find that values of C^'^ out- 
side the interval [—56, 92] MeV fm^ lead to instabilities as 
well. The limits of the stable regions, however, should be 
taken with a grain of salt because their values are sensi- 
tive to the details of the calculation, the mass number of 
the nucleus, or the other parameters of the EDF. In fact, 
even when using an a priori unstable parameterization, 
a finite-size instability might fortuitously remain unde- 
tected in the calculation of a finite nucleus, depending 
on convergence criteria, cutoffs in the numerical repre- 
sentation, the initial conditions of the calculations and 
other numerical choices made. An unambiguous way to 
identify one class of finite-size instabilities is through the 
calculation of the response function of the model system 



of isotropic homogeneous infinite nuclear matter (INM) 
to perturbations of the density in random phase approx- 
imation (RPA) [M]. When an instability occurs at infi- 
nite wavelength, the entire bulk of homogeneous nuclear 
matter undergoes a transition into a different homoge- 
neous state of nuclear matter. This kind of instability 
can be identified from the values of the Landau parame- 
ters discussed in Sect. 



HE 



14, 70, 84J. If instead the in- 
stability occurs at a finite wavelength, the homogeneous 
nuclear matter can undergo a phase transition into an 
inhomogeneous phase, i.e. it exhibits a finite-size insta- 
bility. The former instabilities are driven by the bulk 
terms in the EDF, whereas the latter are driven by the 
terms that contain a nabla or Laplacian acting on a den- 
sity and that are zero in homogeneous INM. Recently, 
the calculation of the linear response of the full Skyrme 
EDF ( 13a|13b ) with tensor terms in force-based [85] and 
general functional [86j frameworks has become available. 
There might also be, however, a second class of insta- 
bilities which are related to surface modes [57]. A more 
detailed analysis of the finite-size instabilities will be re- 
ported elsewhere [55] . 

In a force-base framework and for the forces consid- 
ered here, there is at least one of the C^^ St ■ Asj and 
CY'^ (V • St)^ terms causing non-convergence (see also 
Fig. |3| . To suppress this unphysical behavior, these four 
coupling constants are set to zero in all calculations re- 
ported below unless stated otherwise. 



C. General features 

Before discussing the rotational properties of ^^'^Hg, we 
analyze the evolution of the total energy and the single- 
particle spectra as a function of deformation. 

The total contribution of the tensor terms relative to 
their value at spherical shape is plotted in the panel (a) of 
Fig. [5] For T22, the time-even tensor Ejj contribution is 
close to zero for all deformations. For the other param- 
eterizations, it increases relative at small deformations 
and then follows an almost 'oscillatory' pattern. As can 
be seen from panel (b), the differences between the pa- 
rameterizations seen in panel (a) are strongly attenuated 
in the energy curves. As discussed in Article II, this last 
result comes from an intricate compensation between all 
energy contributions to the EDF. The location and depth 
of the superdeformed minimum around /32 = 0.65 that is 
the key point for the subsequent discussion is very similar 
for all parameterizations. 

The proton and neutron Nilsson diagrams are pre- 
sented in Fig. [6] We recall the conclusion of Article II 
on deformed nuclei that for parameterizations with differ- 
ent strength of the tensor terms the differences between 
the single-particle spectra at sphericity are almost com- 
pensated at large deformation by the changes in slope 
for of the single-particle energies in the Nilsson diagram. 
Hence, the single-particle spectra around the Fermi en- 
ergy for strongly deformed nuclei are often found to be 
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FIG. 5: (color online) Change of the (a) contribution from 
the tensor terms relative to the values at the spherical shape 
and (b) deformation energy relative to the spherical shape for 
^^■'Hg obtained with the parameterizations T22, T26, T44, 
and T62. The energy scale is the same for the two panels. 



very close for all T/J tensor parameterizations in spite 
of significant changes at sphericity. The same property 
is observed in the case of ^^'^Hg. At the Z = 80 superde- 
formed shell gap, the single-particle spectra for all pa- 
rameterizations except T62 lie on top of each other. The 
neutron single-particle spectra follow each other closely 
for T22 and T44 with some small deviations for T26 and 
T62. 



1. SLy4, T22, and T44 

We first compare results obtained with the SLy4, T22 
and T44 parameterizations. 

The dynamical moments of inertia J'^^^ave plotted as 
a function of hio in panel (a) of Figure [t] and the charge 
quadrupole moments Qc in panel (b). The differences 
between the J''^^'> calculated with SLy4 and T22 are 
marginal, even though they correspond to very different 
coupling constants Cf^ and Cf. The J'^^'' obtained with 
T44 increases slightly faster with a plateau appearing for 
a smaller value of huj but, overall, the moments of inertia 
obtained with three interactions present the same behav- 
ior. As discussed in Article II, the position of deformed 
minima may depend on the parameterization but this is 
not the case here as can be checked from panel (b) of 
Fig. [7] and also Fig. [5] The charge quadrupole moments 
Qc values obtained with the three parameterizations dif- 
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02 

FIG. 6: (color online) Proton and neutron Nilsson diagrams 
in ^^''Hg for the T22, T26, T44, and T62 parameterization. 

fer only by about 1% at all spins. This indicates that the 
differences between the moments of inertia are mainly 
due to the differences in the relative weight of the con- 
tributions in the EDF, and not to a change in the shape 
of the nucleus. 

In Figs.[8]and[9) we present the contributions of various 
time-even and time-odd terms to the total energy as a 
function of hut. Their labels refer to their density content 
in the EDF. The total energy is the sum of the kinetic 
{Ekin + Ec.m.), the pairing {Epatr + Eln), the time-even 
{Esk (time even)) and the time-odd {Egk (time odd)) 
Skyrme parts of the EDF. 

The excitation energy at ( J^) = 54 h, which is the value 
at which we stopped the calculations, is 14.1 MeV for 
SLy4 and T22 and 13.6 MeV for T44. AU parts of the 
EDF contribute to this excitation energy in a very simi- 
lar way for all parameterizations, around 7 MeV for the 
pairing energy, 5 MeV for kinetic energy, 2 MeV for the 
Skyrme EDF and 0.5 MeV for Coulomb. The Skyrme 
contribution results from a cancelation between the time- 
even (5 MeV) and time-odd (-3 MeV) parts. 

As expected from the Coriolis-anti pairing effect |52) . 
the pairing energy shows the largest variation with Huj, 
decreasing to less than half its value at spin zero. The 
energies of the kinetic, Coulomb and time-even terms of 
the Skyrme EDF change only by a small fraction of their 
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FIG. 7: (color online) (a): Proton (tt), neutron (i/), and 
total (t) dynamical moments of inertia as a function of the 
rotational frequency for the SD band in ^^^Hg with the SLy4, 
T22, and T44 parameterization, (b): The charge quadrupole 
moment in as a function of the rotational frequency for the 
SLy4, the T22, and the T44 parameterization. 



absolute values when going from (J^) — Oh to (J^) — 
54 h. The time-odd terms of the EDF start from zero at 
spin zero and increase in absolute value with hiu. The 
sum of all time-odd contributions is negative and cancels 
more than half of the contribution brought by the time- 
even terms in the Skyrme EDF. 

Let us first discuss the energy contributions of the 
time-even terms, plotted in Fig. |8] As discussed in Arti- 
cle II, the relative contributions of all terms in the EDF 
differ between each of the TIJ parameterizations. As 
a result, the magnitude of the various energy contribu- 
tions is slightly different for the three parameterizations. 
Their Huj dependence, however, is very similar. By con- 
struction, the tensor contribution Ejj is exactly zero for 
SLy4 and almost zero for T22 at all angular momenta. 
One can note that the energy contribution i?pvj is more 
attractive for T44 than for T22 and SLy4. As discussed 
in Article I, this results from the strong correlation be- 
tween the spin-orbit and tensor coupling constants. In 
fact, the ratio of Epyj{T22)/Ep^j{TM) is very close to 
the ratio of the coupling constants CY-\T22) /C7"'(T44). 

The time-odd contributions to the EDF are plotted in 
Fig. [9) From top to bottom and left to right are shown 
(a) the spin terms Egg, (b) the time-odd terms coupled by 
Galilean invariance to the time-even effective mass i?o,-. 
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FIG. 8: (color online) Evolution of the time-even terms in the 
EDF as a function of the rotational frequency huj for SLy4, 
T22, and T44 for the ground-state superdeformed band of 
"*Hg. 
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FIG. 9: (color online) Same as Fig. [s] but for the time-odd 
terms in the Skyrme EDF. 
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FIG. 10: (color online) Dynamical moments of inertia of 
different terms in the EDF as a function of the rotational 
frequency for SLy4, T22, and T44 in the calculation of the 
ground-state superdeformed band of ^^^Hg. 
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FIG. 11: (color online) Same as Fig. 10 but for the time-odd 
terms in the Skyrme EDF. 



(c) the spin-current tensor E^t + E^p and (d) the time- 
odd spin-orbit Es\j-xj. The sum of the spin terms that 
contribute to the equation of state of infinite homoge- 
neous spin-polarized nuclear matter is shown in panel (e) , 
and the sum of all time-odd contributions in panel (f). 
We recall that the iJ^As and SvsVs terms not shown in 
the figure are set to zero as they cause finite-size instabil- 
ities, see Sect. |III B[ All time-odd terms start from zero 
and change rapidly with increasing hio. The E^t + E^p 
term is zero by construction for SLy4 and is negligible 
for T22. For T44, it decreases down to -800 keV at 
(J^) =54/1. 

All other time-odd contributions are very similar for 
SLy4 and T22. This is not surprising since the coupling 
constants of these terms are very similar. It indicates also 
that the additional Cj'sfTt and Cfs-Ft terms of T22 do 
not introduce a large polarization. The situation is differ- 
ent for T44, for which all time-odd terms but — Cj^jj take 
very different values. These changes are due to the larger 
and Cf coupling constants and the increased spin- 
polarization they induce. All time-odd terms containing 
the spin density, however, tend to cancel each other for 
all parameterizations, such that the sum of all time-odd 
terms is very close to Ejj. 



Using Eq. ( 29 ) , the dynamical moment of inertia can 
be decomposed into various contributions to the EDF 
that are plotted in Figs. 10 and 11 The offset in total 



energy between the parameterizations has no effect on 
the moments of inertia. 

The evolution with Hlj of the energy contributions be- 
ing similar for all parameterizations, the general pat- 
tern of the contributions to J'^^-' is the same for most 
terms. At low huj, the main contribution to J^^^^ is pro- 
vided by the pairing energy Epair + Eln and represents 
about 75%. The Skyrme contribution brings the remain- 
ing 25%, the contributions from Ekin + Ec.m. and Ecoui 
being very small. These two contributions grow rapidly 
with huj reaching about 65% at high spin whereas the 
pairing contributions shrinks to approximately 25%. The 
Coulomb contribution grows slowly with spin but never 
exceeds 10% of the J^^) _ 

The time-even and time-odd contributions coming 
from the Skyrme EDF have opposite signs and to large 
extent cancel each other at high spin as the total con- 
tribution of the Skyrme EDF to J''^^-' does not exceed 
±10%. This cancelation however, is not a generic feature, 
as will be seen with the examples of other parameteri- 
zations and other nuclei, although the £^5^ (time-even) 
and Esk (time-odd) counteract each other in all cases 
encountered. For all parameterizations Ejj is the largest 
time-odd contribution. This was already found in an 
earlier study of the time-odd components of the Skyrme 
EDF l^nj , where it was concluded that the cranking term 
mainly induces a nonzero flow of nuclear matter as mea- 
sured by the Ejj term. 

Concentrating on the tensor contributions, the time- 
even Ejj for T44 is small and varies between -4% at 
low spin and 2% of the J^^) ^t high spin. On a similar 
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scale, the time-odd Est + Esf contribution lowers the 
J^'^^hy at most 7%. Moreover, Ess and Est + Esf tend 
to counteract each other for T44, the total contribution 
Ess + EsT + EsF being less than 1%. As already argued 
in Articles I and II, the presence of the tensor terms has 
an impact through rearrangements of other coupling con- 
stants in the fit and through self-consistency of the HFB. 
This combined effect modifies all time-even and time- 
odd Skyrme J^^"^^ contributions, whereas the pairing and 
Coulomb contributions do not change. 

In Figs. [12] and [13) we present the proton and neu- 
tron quasiparticlc Routhians, i.e. the eigenvalues of Eq. 
(figl), as a function of fiw for SLy4, T22 and T44. All 



Routhians are characterized by their parity and signa- 
ture, by the j-component of the dominant single-particle 
state in the quasi-particle wave function along the axis 
of largest elongation at hw=Q, and by their particle (p) 
or hole (h) character. As discussed before, the proton 
and neutron Nilsson diagrams of T22 and T44 display 

This is reflected in the 
At Sa;= 



only small differences (Fig. l6) 



12 and 13 



=0, the 



quasi-particle spectra in Fig, 
ordering of the lowest quasi-particle states is the same 
for SLy4, T22, and T44, although the values of the en- 
ergy are parameterization-dependent. For the protons, 
the three lowest quasi-particle Routhians are the p 5/2+ 
and the h l/2~ and h 1/2+ . The first significant differ- 
ence between T44 and the two other parameterizations 
is the position of the h 3/2+ state, which is much closer 
to the h 1/2+ level for T44. This position affects the 
alignment of the p 5/2+ Routhian with increasing huj^ 
making it slightly faster for T44 because of a stronger 
interaction with the higher-lying positive-parity states. 
In the same way, for the neutrons, the exchange in po- 
sition of the h 3/2~ and the p 5/2~ quasiparticle states 
at huj—Q has an impact on the alignment of the h 5/2~ 
state. Indeed, the h 3/2~ quasiparticle Routhian occurs 
lower in energy for T44. Since it aligns more rapidly 
than the p 5/2~ state, the interaction with the h 5/2^ 
state takes place at a lower frequency, which affects the 
alignment of this last quasiparticle Routhian. 



2. T26, T44, and T62 

We now proceed to a comparison between results ob- 
tained using the T26, T44, and T62 parameterizations. 
They differ by their value of the isovector C( coupling 
constant while they share the isoscalar Cq = 120 MeV 
fm^ constant. This comparison probes the direct impact 
of a variation of C/, but also its indirect impact due 
to the changes of the other coupling constants that are 
readjusted for each parameterization. 

The dynamical moments of inertia ^j'-^-' and the charge 
quadrupole moments are presented in respectively panel 
(a) and (b) of Fig. 14 respectively. 

Whereas the J^^^btained with the T22 and T44 pa- 
rameterizations (Fig. [7| have a different slope, the be- 
havior of the fQj. rp26^ T44 and T62 is similar up 
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FIG. 12: Proton quasi-particle Routiiians for tiie ground- 
state superdeformed band of ^^''Hg for the SLy4, the T22, 
and the T44 parameterization. The (parity,signature) com- 
binations are indicated in the figure. All Routhians are ad- 
ditionally characterized by the j-component of the dominant 
single-particle state in the quasi-particle wave function along 
the axis of largest elongation at hjjj=0 and by their particle 
(p) or hole (h) character. 
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FIG. 13: Same as Fig. [12] but for neutrons. 



to 0.3/iaj and differs at large spins by the height of the 
plateau. The proton moment of inertia displays a peak 
for T62, which is absent for the other parameterizations. 
This difference can be related to the proton quasiparticle 
Routhians, presented in Fig. 
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In contrast to the other 
parameterizations, it is not the p 5/2+ but the h 3/2+ 
Routhian that is the lowest quasiparticle state for T62. 
The p 1/2+, is at 1.8 MeV for T26 and higher for aU other 
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FIG. 14: (color online) (a) Proton (vr), neutron (v), and total 
(t) dynamical moments of inertia in ^^''Hg as a function of 
the rotational frequency for the SD band in ^^^Hg with the 
T26, the T44, and the T62 parameterization, (b) The charge 
quadrupole moment in ^^^Hg as a function of the rotational 
frequency for the T26, the T44, and the T62 parameteriza- 
tion. 



parameterizations but T62 for which it is only at 1.6 
MeV. With increasing hu, the energy of its negative sig- 
nature partner decreases quickly up to a bending around 
0.4 huj, which is the frequency at which the moment of 
inertia displays a peak for T62. The charge quadrupole 
moments differ by less than 1%, 
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b). 



see Fig. 

The contributions to the dynamical moments of iner- 
tia J''^^' of respectively the time-even and the time-odd 



terms in the Skyrme EDF are presented in Figs. 17 and 
All contributions to the J'*-^-' have a similar behavior 
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-E„ 



for the three parameterizations, except for the Epp - 
terms. 

As for the total moment of inertia, the peaks obtained 
in most terms at large hiu are the most pronounced for 
T62 and the smoothest for T26. Similar to the decompo- 
sition of the J^^^) foj. sLy4, T22, and T44, the time-even 
and time-odd Skyrme contributions are of equal impor- 
tance because they are of the same order of magnitude. 
The more pronouncedly peaked the behavior of the other 
contributions to the J^^^^ the more the Skyrme contri- 
bution tends to counteract it. Indeed, whereas the dif- 
ference in height of the peak between the T26 and the 
T62 parameterization is approximately 60 MeV~^ for 
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FIG. 15: Same as Fig. [12] but for the T26, T44, and T62 
parameterization. 
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contribution, a difference that is even 



FIG. 16: Same as Fig. [15) but for neutrons. 



more intensified by the contribution of Epair + -Eljv and 
Ecoui, the difference between the total J^^^ of the re- 
spective parameterization in Fig. [14] is only about 30 h 
MeV~^. The Skyrme contributions are responsible for 
the decreased difference in peak-height of the J'^^^ . 



D. Time-odd terms 

The parameters of the EDF are usually adjusted to 
binding energies and r.m.s. radii of doubly-magic nu- 
clei and on properties of infinite nuclear matter |47j . As 
discussed in Sect. II B 2 for a force-generated EDF, 
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FIG. 17: (color online) Same as Fig. [lO|but for the T26, the 
T44, and the T62 parameterization. 
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FIG. 18: (color online) Same as Fig. |10[ but for the time-odd 
terms in the Skyrme energy density functional with the T26, 
the T44, and the T62 parameterization. 
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the coupling constants of the time-odd part are en- 
tirely fixed by those of the time-even part although 
they are rarely directly constrained by observables. By 
contrast, the coupling constants of the time-odd terms 
that are not constrained through Galilean invariance 
are a 'priori undetermined in a more general EDF. 
Their adjustment has been the subject of several studies 
26, 28, 31, 32, 40, 55, 89,,^. The possibility of using 
band terminating states to constrain the time-odd terms 
is discussed in Refs. [311 ES] and in Ref. [40] the ef- 
fect of the spin-isopin coupling constants of the Skyrme 
EDFs on predictions for Gamow- Teller distributions is 
investigated. In the latter study, a local fit of the Cf [0], 
Cf [pnm], C^" , and coupling constants to the existing 
data was performed. 

To study the effect of time-odd tensor terms in rotating 
nuclei, we proceed in the following way. Until now, the 
coupling constants of the Sf • Asi and (V • St)^ terms have 
been set to zero for all T/J parameterizations because 
these terms can induce finite-size instabilities. In Sect. 



Ill B we have determined empirical limits between which 
these instabilities do not appear. To maximize the effect 
of these time-odd terms on j'^'^^ , we have taken values for 
C^^ and Cj^* close to their respective limits of stability. 

The other time-odd tensor terms SfTt and SfFj are re- 
lated to time-even tensor terms through Galilean invari- 
ance. It is therefore not desirable to vary their coupling 
constants directly and we have proceeded in an indirect 
way. We have seen in Sect. |III C 1| that the contribu- 
tions of Ess and Est + Egp to J^'^^ act in opposite ways. 
We have therefore changed the values of Cl\pmn\ by ±50 
MeV fm^ around their Skyrme force values at the satura- 
tion density of nuclear matter, which are CQ[p„m] k, 150 
MeV fm^ and C{\pnm\ ~ 100 MeV fm^ for all parameteri- 
zations. In addition, all CI have been taken independent 
of density, choosing its value at the saturation density 
of nuclear matter for the coupling constant that is not 
varied. All variations stay within the stability limits of 
the conditions outlined in [T3] . The discussion of Fig. 26 
below will indeed confirm that a variation of the Cf [poj 
coupling constant has an indirect effect on the St • Tt and 
St • Fi tensor terms. 



1. SLy4 

Let us first limit ourselves to the SLy4 parametrization 
and determine whether the effect of the variation of the 
time-odd coupling constants can be differentiated from a 
pairing effect. The dynamical moments of inertia deter- 
mined with reduced and increased pairing strengths are 
plotted in Fig. [TOj As expected at low spins, a reduction 
of pairing increases i/*-^-* , whereas an increase lowers it. 
The height of the plateau at high spins moves in opposite 
direction but it appears for similar values of huj. 

In Fig. 20 we show the dependence of J''-^-' on the 
isoscalar Cn and isovector Cf coupling constants. The 



only noticeable change in JZ^^^ is a slight shift of the 
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FIG. 19: (color online) Dynamical moments of inertia as a 
function of the rotational frequency for different choices of the 
strength V — Vp = V„ of the DDDI pairing interaction, Eq. 
(18 1. The SLy4 parameterization was chosen in the particle- 
hole channel of the interaction and the pairing strength is 
given in units MeV fm~^. 
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FIG. 20: (color online) Dynamical moments of inertia as 
a function of the rotational frequency for different density- 
independent choices of the C" coupling constant. If the 
isoscalar Cq is varied, then the density-independent Cf is cho- 
sen to be equal to Cf [nm] and vice versa. All other coupling 
constants in £sk are determined by the SLy4 parameteriza- 
tion and the Ct coupling constant is given in units MeV fm^ 



plateau for Cq. The effect is small but clearly different 
from the effect of a variation of the pairing strength; also, 
it is larger for t = than for < = 1. 



2. T22, T26, TU, T62 

We n ow tur n to interactions including a tensor term. 
In Figs. 2l|24 we present the variations of J''^^^ as a func- 
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TABLE I: Spin and spin-isospin Landau parameters for SLy4, 
T22, T26, T44, and T62 parameterizations. 
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0.904 
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-0.044 


0.870 
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0.526 
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0.415 


-0.253 
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tion of (a) Cf ^ (b) C^"", and (c) for T22, T26, T44, 



TABLE II: Spin and spin-isospin Landau parameters go and 
g'o for the density independent variations of the Ct [0] coupling 
constants of SLy4, T22, T26, T44, and T62 described in the 
text. The units of the Ct[0] coupling constants is MeV fm'^. 



and T62 respectively. In the panels (a), one can see that 
the dynamical moment of inertia presents a significant 
decrease at high hut for a positive value of Cj^" {t = 
and t = 1) and is nearly unaffected by the Sf • Asj for 
a negative value. The result is inverted for the (V • St)^ 
term, (panels (b)), for which the coupling constant has 
to be negative to have a visible effect. 

This behavior can be understood by looking to Fig. 
|4j which is devoted to C^^ but is representative for all 
four coupling constants. The energy of the Sq ■ Asq term 
varies rapidly as a function of C^" when it is positive 
and close to the value leading to instabilities. Owing 
to self-consistency effects, several other time-odd terms 
vary rapidly for values of C^^ close to its maximal value 
before instabilities appear. For negative values of C^*, 
the corresponding energy behaves much more smoothly 
and therefore does not affect significantly the moment of 
inertia. The picture is the same for C^" and and for C^^, 
except that then, rapid changes of some energy terms are 
obtained for negative values of the coupling constant. 

The panels (c) show how J'^^'> is affected by variations 
of C^. The changes are very similar for T22 and T26 
and much larger for T44 and T62 when the t = cou- 
pling constant is varied. Naively, one would expect that 
the changes of the J^^^' with respect to the "original" 
J'^'^^ (obtained with the force-generated value of C^) 
should be similar for all parameterizations. Our result 
clearly indicates that the dynamical moment of inertia is 
not sensitive to the individual values of the various cou- 
pling constants but rather to some specific combinations 
of them. 

The values of the Landau parameters (see Sect. II E| in 
the spin and the spin-isospin channels are given in Tablcjl] 
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FIG. 21: (color online) Dynamical moments of inertia as a function of the rotational frequency for variations of (a) C^" , (b) 
the CY" , and (c) Ct coupling constants. All coupling constants of the Skyrme EDF not explicitly mentioned in the legend are 
determined by the T22 parameterization. For the density- inde pendent variations of the C| coupling constant (c), the same 
prescription as for SLy4 and explained in the caption of Fig, 
in units MeV fm 
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is followed. The CY" and Ct^" coupling constants are expressed 
and the coupling constant is given in units MeV fm^. 





for all parameterizations considered in this work. Table presents the and Landau parameters corresponding 
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gs. 21 


24 


21|24 


(c) 



to the variations of that we consider in Fi gs. [2T 
A close inspection of these Tables and of F 
puts in evidence some clear trends: 

• For a given interaction, an increase (decrease) of 
go increases (decreases) the value of hu where the 
moment of inertia is maximum. It also broadens 
(sharpens) the peak of the J^^^ curve. 

• Similar values of go lead to a comparable depen- 
dence of J'^^^ on huj, as is e.g. the case for T22 
and T26 for the three values of go given in the Ta- 
bles but also for T22, T26, and T44 (C^ = 200 MeV 
fm3), which have 50 ~ 0.9, or for T44 (C^ = 100 
MeV fm^) and T62 where go ^ -0.1. 

• The behavior of J^^^) 

as a function of hu does not 
depend much on go values larger than 0.4. Its vari- 
ation as of a function of go is much larger for go 
lower than 0.4. 

This confirms and extends a conclusion of Bender et 
al. |40j . These authors have shown in a similar man- 
ner and for another Skyrme parametrization SkO' that 
the dynamical moment of inertia J''^^^ of the SD band 
of ^^^Dy depends mainly on the spin-isospin Landau pa- 
rameter (7q and not so much on the actual values of the 
individual coupling constants Cf and Cf . Different com- 
binations of EDF coupling constants leading to the same 
value of g'o were found to change little to the J^^^ Note 
that these authors have also varied C^'^, extending their 
study well outside the limits for which we find instabili- 
ties. This discrepancy might be due to the use a different 
technique for solving the mean- field equations, i.e. by 
means of an expansion on an oscillator basis instead of 
using a cartesian mesh. 

Note also that empirical values for go ^ 0.4 [21] and 
1-4 < g'o < 1-6 [92H94] have been derived from Ml and 
Gamow- Teller response. The values quoted in Table |l] 
and |lT] are in all cases except one very different from the 
empirical values. 
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FIG. 25: (color online) Contributions to the dynamical mo- 
ment of inertia coming from the time-even terms of the EDF 
as a function of the rotational frequency for variations of the 
CI coupling constant. All coupling constants not explicitly 
mentioned are determined by the T44 parameterization. The 
CI coupling constants are expressed in units MeV fm^ 



The previous analysis shows that the , Cf and Cj ^ 
coupling constants are linked by their contribution to the 
values of the Landau parameters go and g'^. One might 
wonder whether variations of affect mainly the terms 
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FIG. 26: (color online) Same as Fig. 25 but for the time-odd 
terms in the Skyrme EDF. 



IV. RESULTS FOR THE SUPERDEFORMED 
BAND IN i^^DY 

A. General comments 

The superdeformed rotational bands known in ^^^Dy 
exist in a very different regime than the one of ^^^Hg. 
SD band in the A w 150 region have only been detected 
for spins higher than 20S. For such angular momenta, 
pairing correlations are strongly weakened by the Cori- 
olis anti-pairing effect. Hence, pairing is expected to 
play only a minor role in that region. Early studies 
by Bengtsson et al. [SS] using the Nilsson-Strutinsky 
approach have demonstrated that the behavior of the 
SD bands is strongly influenced by the number of nu- 
cleons that occupy the intruder orbitals. SD bands in 
the Dy-region have been studied extensively within self- 
consistent mean-field approaches [511 SOI US] • More- 
over, they have been used as a laboratory to study the 
time-odd terms in the EDF [28 l l40 l [55] . 



B. General features 



in the EDF depending on the spin-density Sj. To ana- 
lyze how all terms in the EDF are affected by a variation 
of C|, we have decomposed the J'^^-' obtained for the 
modified T44 parameterizations in the same way as in 
Figs. 25 and 26 Differences appear only for large values 



of hoj. At hu}=OA, taking = 100 MeV fm^ modifies 
the total energy obtained with the T44 parameterization 
by about 600 keV. This change results from the partial 
cancelation of larger changes with different signs of all 
terms of the EDF, including the pairing and Coulomb en- 
ergy. In particular, the Skyrme time-even and time-odd 
contributions decrease the energy by about 900 keV and 
500 keV respectively. The energy differences obtained for 
Cq = 200 MeV fm^^ are of the same order of magnitude, 
but are much smaller for variations of C(. These energy 
changes at high hcu affect the slope of the different en- 
ergy contributions and therefore the ir con tributions to 



Ulel One clearly 



the ^ xhese are plotted in Figs, 
sees that the modification of the coupling constant of a 
relatively small term, in casu the E^s, affects all other 
terms through self-consistency and how little changes in 
the energy can make a large difference in the J'^^\ 

Our analysis demonstrates that the variations of pair- 
ing and of the time-odd terms have clearly distinguish- 
able effects on the shape of the J'^'^^ . The time-odd terms 
influence the slope of the J^*^^^ and the hu at which the 
plateau occurs. While the effect of variations of the C^'^ 
and Cj^'* coupling constants on the J'*^^^ is rather small 
and depends on the sign of the coupling constant, varia- 
tion of the Cq coupling constant may lead to significant 
changes in the 



The dynamical moments of inertia (panel (a)) and the 
charge quadrupole moments (panel (b)) are shown in Fig. 



27|for SLy4, T22, and T44 and in Fig. [28] for the T26, 
T44, and T62. As for ^^^^Hg, the calculated with 

SLy4 and T22 are very close, with small differences at 
low spin. The main difference between the parameteri- 
zations is the presence of a strong peak at low spin for 
T62, which is less pronounced for T44 and absent for the 
other parameterizations. In all cases, the difference be- 
tween the J'^^^ is caused by the neutrons. The Qc values 
presented in the panels (b) differ by less than 2% and all 
display the same behavior, increasing until huj—0.6 McV 
after which they slowly start decreasing again. 

The total energy and the J^*^^^ are decomposed into 
thei r time- even and time-odd components in Figs. |29|30| 
and 3l|[32 respectively, as it was done for ^^^Hg. Decom- 
positions are only shown for the SLy4, T22 and T44. The 
results for T26 and T62 are very similar. The difference 
of energy between the lowest {{Jz} = 28fi) and highest 
{{Jz) = 74fi) states that we have calculated amounts to 
26.3 MeV for SLy4 and T22 and 25.5 MeV for T44. All 
parts of the EDF give contributions of the same order of 
magnitude to this change in energy, about 12 MeV for 
the Skyrme EDF, 10 MeV for the kinetic energy, 7 MeV 
for the pairing energy and —2 MeV for the Coulomb en- 
ergy. The Skyrme contribution can be decomposed into 
18 MeV from the time-even and —6 MeV from the time- 
odd terms. 

Focussing on the time-even contributions (Fig. 



291, 



their huj dependence is very similar for all parameter- 
izations, which is also reflected in the c orre sponding 

As ex- 



31 



contributions to i/*-^-* presented in Fig. 
pected, the pairing energy is smaller than in the A w 
190 region, because the superdeformed bands in the 
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FIG. 27: (color online) (a) : Proton (tt), neutron (v), and 
total (t) dynamical moments of inertia as a function of the 
rotational frequency for the SD ground band in ^^^Dy with 
the SLy4, the T22, and the T44 parameterization . (b) : 
The charge quadrupole moment in ^^^Dy as a function of 
the rotational frequency for the SLy4, the T22, and the T44 
parameterization. 



FIG. 28: (color online) (a) : Proton (tt), neutron (v), and 
total (t) dynamical moments of inertia as a function of the 
rotational frequency for the SD ground band in ^^^Dy with 
the T26, the T44, and the T62 parameterization, (b) : The 
charge quadrupole moment in ^^^Dy as a function of the ro- 
tational frequency for the T26, the T44, and the T62 param- 
eterization. 



A 150 region only occur at higher spin. The ratio 
Epv J {T 22)/ Epy J (TM) is again approximately equal to 
the ratio of the corresponding coupling constants. For 
the time-odd terms (Fig. 30), the contributions of Est 
and EsF terms that appear for T44 cancel out the Ess 
term, which is much larger without the inclusion of a 
tensor term. 

The time-even and time-odd contributions to the J''^-' 
(Figs. 31p2 ), indicate that the Esk (time-even) is the 
largest contribution at low huj, amounting up to 96% of 
the total value for T22 and about 80% for T44. The 
pairing energy Epair is the second largest contribution, 
averaging about 70%, and is mostly canceled out by Ecoui 
(about -40%) and the time-odd Esk terms (about -25%). 
The Ekin + Ec.m. contribution is negligible at low Hui . 
With increasing huj, the Epair + E^j^i and Ecoui con- 
tributions drop to 13% and 2% respectively, while the 
Ekin + Ec.m. contribution quickly grows to about 50%. 
Whereas the time-odd Egk contribution stays more-or- 
less constant with Hui , the time-even part of Esk reduces 
to about 70%. Hence, in contrast to ^^''Hg, the Skyrme 
EDF plays a more important role in the decompositions 
of the total energy and of the J"'^^ . Even though changes 
the time-even and time-odd components of the Skyrme 



EDF no longer almost cancel out, they still counteract 
each other. 

Finally, the proton and neutron quasiparticle Routhi- 
ans for the SLy4, T22 and T44 parameterization are pre- 
sented in Figs. [TSpB) Again, the low-lying quasiparticle 
Routhians are very similar for all parameterizations and 
subtle differences are observed in the alignment of the 
Routhians. The main difference observed between the 
parameterizations is the location and evolution of the 
neutron h 1/2+ state. When going from SLy4 to T22 
and further on to T44, the h 1/2+ quasiparticle state 
starts of at lower energy and the minimum becomes more 
pronounced. 



V. DISCUSSION AND CONCLUSION 

We have studied the impact of the introduction of ten- 
sor terms in the Skyrme energy density functional on the 
dynamical moments of inertia of superdeformed bands. 
The excellent description of these bands by conventional 
EDFs was a major success of microscopic mean-field 
models in the 90s. Therefore, it is important to verify 
that the inclusion of a zero-range tensor force does not de- 
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FIG. 29: (color online) Evolution of different terms in the 
energy density functional as a function of the rotational fre- 
quency for the SLy4, T22, and T44 parameterizations in the 
calculation of the ground state superdeformed band of ^^^Dy. 
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FIG. 30: (color online) Same as Fig. 29 but for the time-odd 



FIG. 31: (color online) Dynamical moments of inertia of dif- 
ferent terms in the energy density functional as a function of 
the rotational frequency for the SLy4, T22, and T44 param- 
eterizations in the calculation of the ground state superde- 
formed band of ^^^Dy. 
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FIG. 32: (color online) Same as Fig. |31[ but for the time-odd 
terms in the Skyrme energy density functional. 
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FIG. 33: Proton quasi-particle Routhians for the ground- 
state superdeformed band of ^^^Dy for the SLy4, the T22, 
and the T44 parameterization. The (parity, signature) com- 
binations are indicated in the figure. All Routhians are ad- 
ditionally characterized by the j-component of the dominant 
single-particle state in the quasi-particle wave function along 
the axis of largest elongation at hw=0 and by their particle 
(p) or hole (h) character. 
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FIG. 34: Same as Fig. |33| but for neutrons. 



stroy this agreement. In our study, special attention was 
paid to the time-odd tensor terms in the Skyrme EDF, 
which are zero when self-consistent time-reversal symme- 
try is not broken, as is the case for the ground states of 
the spherical and deformed nuclei that were studied in 
Articles I and II. 

In order to disentangle their respective role, we have 
tested a selection of four T/J parameterizations, intro- 



duced in Article I, which represent a wide range of values 
for the isoscalar and isovector tensor coupling constants. 
As a reference without tensor terms, we have included the 
SLy4 parameterization. We have found that the inclu- 
sion of tensor terms in the Skyrme EDF does not change 
the overall behavior of the dynamical moments of inertia, 
although differences in slope and location of plateau do 
occur. This results from an intricate compensation mech- 
anism due to the self-consistency that is implemented at 
two different levels in our method: in the fitting proce- 
dure of the interactions considered and in the solution of 
the mean- field equations. 

Even though the energy contribution of the time-even 
Ejj tensor terms is in general an order of magnitude 
larger than that of the time-odd E^x + E^p tensor terms, 
the time-odd tensor terms evolve more rapidly as a func- 
tion of rotational frequency such that their contribution 
to j/'^-* is of the same order of magnitude as the one of the 
time-even tensor terms. Similarly, the Skyrme time-even 
and time-odd energy contributions Esk typically differ by 
three orders of magnitude but have similar contributions 
to J"^^^. In all cases encountered, the time-odd Esk con- 
tributions to the j/^^' partially cancels out the time-even 
Esk contribution. 

A detailed study of the time-odd terms in the Skyrme 
EDF has shown the following features: 

1. We have seen that the values of the coupling con- 
stants of the time-odd tensor terms that contain 
derivatives of spin densities (£'sAs and i^vsVs) have 
to be chosen within strict limits to avoid finite-size 
instabilities. Such instabilities were encountered 
for all T/J parameterizations considered. There- 
fore, we adopted the functional point of view in 
our study and put their coupling constant to zero. 
By contrast, these instabilities are not encountered 
in spherical QRPA studies using the same TJJ pa- 
rameterizations [TTI - fTT] , presumably because of the 
non-variational character of QRPA. A point of spe- 
cial interest will be the analysis of finite-size insta- 
bilities using the technique of Ref. [SH |85] that is 
currently underway. 

2. The effect of modifications in the strength of the 
pairing interactions can be clearly distinguished 
from modifications of the time-odd terms in the 
Skyrme EDF that are not restricted by symmetry 
considerations. 

3. The energy changes that are introduced through 
the presence of the time-odd tensor terms E^x + 
EsF are partially cancelled out by the presence of 
the other time-odd terms that contain the spin den- 
sity St. Consequently the time-odd Esk is almost 
entirely determined by the Ejj contribution. Varia- 
tions of the coupling constant indicate that only 
combinations of the Egx + Esf and E^s terms with 
similar values of go and/or q'q act independently. 

When comparing the dynamical moments of inertia 
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obtained for an EDF including tensor terms to the ex- 
perimental ones, the agreement for T22 and T26 is as 
satisfactory as for SLy4. The T62 parametrization that 
acts only between neutrons and protons in spherical 
symmetry exhibits a pronouncedly peaked behavior that 
does not compare well with the experimental j'^'^^ . Fi- 
nally, the results obtained with the T44 parameteriza- 
tion, which was one of the two parameterizations lead- 
ing to the best results for low-lying collective states in a 
RPA calculation of ^o^pb and ""^Ca [l7], are found to be 
in good comparison with the experimental results. 
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Appendix A: Coupling constants of the Skyrme 
energy density functional in the isospin and the 
proton-neutron formulation 



In Table [Till we provide the relation between the cou- 
pling constants of the Skyrme EDF in the isospin formu- 
lation ( 13a||13b ) and those appearing in the Skyrme EDF 
in the proton- neutron formulation (21). 



Appendix B: Densities and currents in cr8 

In this appendix, we provide the expressions of the 
densities and currents as they are implemented in our 
cranked Hartree-Fock-Bogoliubov solver cr8. This ex- 
tends the discussion of Ref . [37j by the densities and cur- 
rents entering the tensor terms. 

The cr8 code uses a coordinate-space representation 
of the wave functions and fields. The HFB equations are 
solved with the so-called two-basis method, where in an 



iterative scheme the HFB Hamiltonian ( 19 ) is diagonal- 
ized in a single-particle basis that converges towards the 
eigenstates of the mean-field Hamiltonian /i, Eq. (20 1. 



The densities needed to construct the local fields are cal- 
culated in the canonical single-particle basis, which is ob- 
tained by diagonalization of the density matrix p. For a 
detailed discussion of our method of solving the cranked 
HFB equations we refer to Refs. 1^ . 

The cr8 code assumes triaxial symmetry of the nu- 
cleus, where all single-particle wave functions have a 
plane reflection symmetry about the x = 0, y = and 
z = planes. There are several possible choices to 



TABLE III: Coupling constants in the isospin formulation of 
the EDF as a function of the coupling constants in the proton- 
neutron formulation of the EDF in the format C = o-i^i-, 
where the factors a,; are given in the Table. 
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achieve this [97 . The cr8 code chooses the single-particle 
wave functions <I>fc(r,(T) to be eigenstates of 

(i) parity 

P$fe(r, ct) = ^k{-r, ct) = pk^k{r, ct), Pk ^ ±1 : 



(Bl) 



(B2) 



(ii) z signature 

i?,$fc(r,CT) =e"-^"=$fc(r,CT) , 

= M/fe$(r,CT), 77fe=±l, 

(iii) y T-simplex 



Sj$fe(r,CT) = fPi?,$fc(r,CT) , 

= $fc(r,CT), (B3) 

where T is the time-reversal operator. 

A wave function is completely determined by four real 
functions '^k,a (a = 1, . . . , 4) that correspond to the real 
(Re) and imaginary (Im) parts of the spin-up and spin- 
down components of A different numbering of these 
four components was adopted for wave functions of pos- 
itive and negative signature in [27) 
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TABLE IV: Parities of the components "i/k.a of a wave func- 
tion of parity pk with respect to the x — 0, y — 0, and 
z = planes. 
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(B4) 



-1. (B5) 



This choice ensures that each of the four real functions 
^k,a has the same definite reflection symmetry about the 



X, y and z planes, listed in Table IV independently of its 
signature. 



In our code, the local densities and currents (2a|-(|2f[) 
entering the Skyrme EDF are constructed in the canon- 
ical basis. There, they can be expressed as 



p,(r)=Ei;2$t(r,a)$fc(r,a) , (B6a) 
r,(r) = E [V$fc(r, • V$,(r, a) , (B6b) 

c) 



k.(7 



k.a 



-[V1>fc(r,a)]t$fe(r,a)} , 



(B6d) 
3e) 



T,(r) = E [V<i>fc(r,a)]t . [V$fc(r,a')] 



(B6f) 



+ [V$fc(r,a)]t[V-<,'$fc(r,a')]} , (B6g) 

where i;^ are the occupation probabilities and fi. u — 
x,y,z. Expressed in terms of the functions the 



scalar local densities take the form 

4 



k a=l 



(B7a) 
(B7b) 



k a=l 

whereas the vector densities are given by 
j(r) = E "fe(*fc,i V*fc,2 - */c,2V1'fc,i 





= E2t'fe(*/c,l«'fc,3 + 1'fc,2*fe,4) , 
k 


(B8a) 
(B8b) 




= E2t'fe??(^'/c,l*fc,4 - *fc,2^'fe,3) , 
k 


(B8c) 


Sz{r) 


= E '^fe^(*M + *M - *M - *m) , 
k 


(BSd) 




^Y.^vl{[V^k^]-['^^k^3] 

k 






= E2^fc^fc{['^*M]-[V«'M] 

fc 

= J2^hk{[y^kA]-mk,i] 

k 


(B8e) 




(B8f) 



-[V^'fc,4]-[V^'fc,4]}, (B8g) 

- E"fe {0fc,l[V,.*fc,l] + 0fc,2[V^*fc,2] 

fc 

+ '^fe,3[V^^'fc,3] + 0fc,4[V^*fc,4]} , (B8h) 

where we defined the spinor 



/ (t}k,l\ ( V:,*fc,3 + vi^y-^kA + V^*fc,l)\ 

(t>k,2 ^ V:r*fc,4 - ?7(Va*fe,3 - V^*fe,2) 
(Afc,3 V:r*fc,l - ?7(Va*fe,2 + V^^'fe,3) 

V<^fe.4/ VV:,*fc,2 + ?7(Va*fe.l - V^*fc,4)/ 

(B9) 

Finally, the spin-current tensor densities are 

k 

+ «'fe,3[V^*fe,2] - *fe,4[V^«'fca]} , (BlOa) 

Jf^V =E"fe - *fc,l[V^^'fc,3] - 1'fe.2[V^«'fe,4] 

fc 

+ «'fc,3[V^*fc4] H- *fe,4[V^«'fe,2]} , (BlOb) 
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TABLE V: Parities of the nucleon densities with respect to 
the X — 0, y — 0, and z — planes 
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depends also on the angle between the direction q = 
q/|q| of the momentum transfer q = qi — q2 and the di- 
rection of the particles' spins. Because the single-particle 
momenta are restricted to the Fermi surface, the param- 
eters F' , G, G' , H, H' depend only on the angle 
between qi and q2 and can be expanded into Legendre 
polynomials, i.e. 

^-^/£^^(cos0), (C3) 
and similar for the other parameters. 



2. The Landau parameters derived from the 
Skyrme energy density functional 



- ^ 



k,3 



fe.4j 



]}. (BlOc) 



for fi = X, y, z. The symmetries (Bl), (B2| and (B3) of 



the single-particle wave functions impose reflection sym- 
metries on the components of the local densities and cur- 
rents, which are listed in Table fVl 



Appendix C: The Landau-Migdal interaction 

1. General considerations 

Landau theory for normal Fermi liquids provides a 
framework for the study of the long-wavelength response 
of a many-body system jj72j |9B]. In this framework, 
the residual interaction is provided by the so-called Lan- 
dau interaction. It determines the response of the sys- 
tem but cannot be used to calculate its ground state. 
Based on this, Migdal developed the Landau-Migdal the- 
ory of finite Fermi systems and applied it successfully to 
study collective modes in atomic nuclei [5S|. The resid- 
ual Landau-Migdal interaction acts between two par- 
ticles with momenta qi and q2 at the Fermi surface, 
|qi| = |q2| = ^F, and is given by 

Wres (qi,q2) = 

7Vo{i^(qi,q2) + -F'(qi,q2) (ti • t^) 

+ G'(qi, q2) (cTi • 0-2) + G"(qi, q2) (cti • (T2)(ti • T2) 



[i?(qi, q2) + i/'(qi,q2) (xi • T2)] 5i2(q)} 



(CI) 



The normalization factor is defined as the average level 
density A*o = Ikpvn* jf^'n'^ at at the Fermi momentum 
kp = (Itt^po)^''^: with ttIq being the isoscalar effective 
mass. The tensor operator 

5i2(q) = 3 (o-i • q) (cr2 • q) - • ^73 (C2) 



To establish the relationship between the Landau pa- 



rameters in (CI) and the coupling constants of the 



Skyrme EDF in symmetric nuclear matter, we follow the 
procedure outlined in [7TJ [SI] . Starting from the Skyrme 
energy density functional, the residual particle-hole in- 
teraction in coordinate space is obtained as 



(riCTiT{,r2Cr2T2l^ph|l"lCriTi,r2(T2T2) = 



(5p(ricriri,r'icrir{)(5p(r2(T2r2,4cr2r^) ' 



(C4) 



From this, the momentum-space matrix elements of the 



residual Landau interaction (CI I are obtained through 
the substitutions Vi = iqi, V2 = zq2, = — «qi, and 
V2 — — iq2 in the Landau limit where initial and final 
momenta qi and q2 are both on the Fermi surface. This 
results in momentum-space matrix elements of the form 



(qi,q2|V'pft|qi,q2> = 

{w^"(qi,q2) + VF^''(qi,q2)(xi-T2 



E^/^-'(qi'q2)(^l-^2)^!V2^}. (C5) 



with 



W/-(qi,q2)=2Co^[po]+4'''^" 



dpi '° 



+ Q(qi-q2)' 
VF''^(qi,q2) =2Cf[po] +Cr(qi -qa)' 

M^-(qi, q2) = + Co^(qi - q2)')<5^. 

+ C'^Cqi - q2)M(qi - 



(C6) 



(C7) 



W^;."(qi,q2) = (2Ci^[po] + Ci^(qi 



- q2j 



+ (qi - q2)^(qi - qz). ■ (C8) 
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The tensor operator 5*12(9) is easily recognized in (C5| 



when we rewrite WJ^^ and 



w^;:(qi,q2) = 




1 ^ P \ / 

-f 3C^o )(qi-q2 










(C9) 


W^™(qi,q2) = 




+ ^C^f)(qi-q2 










(CIO) 



Because qi and q2 arc both on the Fermi surface, 
(qi ^ q2)^ can be rewritten as 2kF{l — cosd), with 9 
being the angle between qi and q2. A straightforward 
comparison between (|C1[) and ( C5 ) then finally gives us 



the relation between the Landau parameters and the cou- 
pling constants in the Skyrme EDF 
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